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Imaging the Molecular Disk Orbiting the Twin Young Suns of 

V4046 Sgr 

David R. Rodriguez^, Joel H. Kastner^, David Wilner^, & Chunhua Qi^ 

ABSTRACT 



We have imaged the disk surrounding the nearby {D ~ 73 pc), ~12 Myr, 

classical T Tauri binary system V4046 Sgr with the Submillimeter Array (SMA) 

CN ! at an angular resolution of ~2". We detect a rotating disk in ^^C0(2-l) and 

' ^^C0(2-l) emission, and resolve the continuum emission at 1.3 mm. We infer 

■^ ^ disk gas and dust masses of ~110 and ~40 Earth masses, respectively. Fits to 

^y\ [ a power-law disk model indicate that the molecular disk extends to ~370 AU 

r^ I and is viewed at an inclination of between ~33° and ~39° for dynamical stellar 

Ph| masses ranging from 1.8 Mq down to 1.5 Mq (the range of total mass previously 

O I determined for the central, 2.4 day spectroscopic binary). This range of disk 

^ I inclination is consistent with that assumed in deducing the central binary mass 

(i.e., 35°), suggesting that the V4046 Sgr binary system and its circumbinary, 

molecular disk are coplanar. In light of the system's age and binarity, the presence 

>• ' of an extensive molecular disk orbiting V4046 Sgr provides constraints on the 

Qv ' timescales of processes related to Jovian planet formation, and demonstrates that 

0^ ■ circumbinary Jovian planets potentially could form around close binary systems. 

(^ ■ Subject headings: circumstellar matter — stars: binaries, individual (V4046 Sgr) 

O . — stars: pre-main sequence — stars: planetary systems — radio lines : stars 

> 

X 

^ • 1. Introduction 

Circumstellar disks around young stars serve both as the sources of material for accreting 
young stars and as the sites of nascent planets orbiting such stars. Numerous studies of such 
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disks have exploited photometry and spectroscopy of dust emission, as manifest in the form of 
a thermal infrared excess above the stellar continuum, to ascertain fundamental properties 
such as disk diraensions, density structu re, dust mass, and dust grain composition (e.g., 
ZuckermanI 1200 ll : iDuUemond et al.l 120071 and references therein) . Observations that can 
establish the composition and evolution of the gaseous component within circumstellar disks 
around pre-MS and young MS stars, meanwhile, are essential if we are to understand the 
processe s involved in Jovian planet formation and the origins o f comets and Kuiper Belt 
objec ts (ILuu fc JewittI l2002l : iHubbard. Burrows, fc Lunind l2002l : lEhrenfreund fc Charnley 
20001 ) . For example, the near-IR lines of rovibrational H2 and CO generally probe the 
terrestrial planet formation zones ( ~2 AU) of gas disk s , interior to the regions where Jovian 



planets are e xpected to form (e. g., ISalyk et al.l 120071 : iNajita et al.l 12008 



20081 : iNaiitaet al. 



howeyer. see iBrittain et al.ll2009l ). Mid-IR pure rotational lines of HoJ 



2OIOI ) and rovibrational lines of H2O, C2N2, and OH (JGlassgold et al. 



Barv et al. 



Bitner et al. 



2008: 



2007, 



2009|, 



and references therein) are potentially diagnostic of molecular gas mass and excitation in 
the innermost Jovian planet-forming zones, and also provide limited information concerning 
the structure, dynamics, and chemistry of these regions. 

Sensitive measurements of (sub)millimeter-wave radio emission from orbiting CO (after 
H2, the most abundant molecular species) around "isolated," young (age ~5 -20 Myr), nearby 
(D 5, 100 pc) stars — with TW Hya as ar chetype — are key in this regard (jZuckerman et al. 



I995I : iKastner et al.lll997l : iThi et al.ll2004j ) . Such stars are old enough th at Jovian planets may 



already be forming or have formed in their circumstellar disks (e.g., iPapaloizou fc Nelson 



2OO5I : iThommes et al.l |2008| ) , and they are close enough for intensive follow-up study via 



interferometric molecular line imaging surveys. Though CO detections of the disks around 
pre-main sequence stars of such advanced ages are exceedingly rare thus far, (sub)millimeter- 
wave molecular imaging studies of the handful of known ~10 Myr-old molecular disks that 
lie within ~100 pc of Earth have yielded a wealth of information concerning the structure 



and e volution of gaseous, planet-forming disks (e.g.. lOi et al.ll2004 l2006l . I2OO8I : iHughes et al. 
2008h . 



V4 046 Sgr is one su ch isolated, "old," nearby star system (age ~12 Myr, distance 
~73 pc; iTorres et al.ll2008l ). The V4046 Sgr system is especially noteworthy in that it consists 
of a close (~9 Rq separatio n, 2.4-day period) pair of ~0.9 Mq stars that are evidently 
still actively accreting mass ( Stempels fc Gahmll2004j . and references therein). A previous 
molecular emission line survey with the Institut de Radio Astronomic Millimetrique (IRAM) 
30-m telescope establi shed the presence of C O, HCN, CN, and HCO^ within the dusty disk 
encircling the system (IKastner et al.ll2008bl ). This makes V4046 Sgr only th e fourth known 
example (after TW Hya, 49 Cet, and HD 141569; IZuckerman et al.l Il995l ) of a pre-main 
sequence star within ~100 pc that possesses a molecular disk; V4046 Sgr is the only known 
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close-separation binary among the m (HP 141569 has a pa ir of M-type companions, ~800 AU 



from the disk-bearing primary, see lWeinberger et al.ll2000l ). To confirm that the CO emission 



detected toward V4046 Sgr with the IRAM 30-m indeed arises in an orbiting, circumbinary 
disk, and to estabhsh or constrain fundamental system parameters such as disk dimension s 



and inchnation, we observed V4046 Sgr with the SubmiUimeter Array (SMA: lHo et aLll2004j ). 



2. Observations 

We obtained data in the SMA's extended (28-226 m basehnes) and compact (6-70 m 
basehnes) configurations on 23 Feb. and 25 April 2009, respectively. The weather conditions 
were excellent on both nights. During each track, the receivers were tuned to provide simul- 
taneous coverage of the ^^C0(2-l), ^^C0(2-l), and C^^0(2-l) transitions as well as adjacent 
(1.3 mm) continuum, with the spectral backends configured to provide velocity resolution of 
0.26 km s~^. The data were edited and calibrated using the MIR software package). The 
passband response was calibrated using the quasars 3C 273 and 1924-292 for the extended 
configuration and 3C 84, 1924-292, and Uranus for the compact configuration. The antenna 
gain calibration was performed using 1924-292 (extended) and 1733-130 (compact). The flux 
scale was determined by bootstrapping observations of 3C 273 and Uranus and is accurate 
at the ~10% level. The continuum was determined using the line-free channels in both up- 
per and lower sidebands for both configurations. From these data we determine a 1.3 mm 
continuum flux of 360 ± 36 mJy. 

We combined the extended and compact configuration data to produce a synthesized 
beam size of 2.2" x 1.5" at the 230 GHz frequency of ^^00(2-1). Similarly, for ^300(2-1) 
we produced a synthesized beam size of 2.3" x 1.7" at its 220 GHz frequency. No C^^O 
219 GHz emission was detected with a 3a upper limit of 0.3 Jy beam~^ (synthesized beam 
size was 2.4" x 1.7"). The standard tasks (CLEAN algorithm) of deconvolution and image 
restoration were performed with the MIRIAD software package. Images of the continuum and 
^^CO were created using Briggs weighting (robust=0 and 1, respectively) for the visibilities. 
To reduce image noise, the ^'^CO and C^^O images were generated using natural weighting. 
The complete observing parameters are listed in Table [TJ 



^See http://cfa-www.harvard.edu/~cqi/mircook.html 
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3. Results 

3.1. CO line maps 

The resulting sequences of velocity-resolved SMA ^^C0(2-l) and ^^C0(2-l) images are 
displayed in Fig. [1] The contours start at 3 and 2 times the rms noise level, respectively, 
and increase as described in the figure caption. These images vividly demonstrate that 
the CO emission from V4046 Sgr arises from a rotating disk. The radial extent of the 
molecular disk is ~5", or ~370 AU at the estimated 73 pc distance to V4046 Sgr. The 
similar overall extent of the CO emission in the N-S and E-W directions, coupled with the 
profound change in emission morphology with radial velocit y, indicates that the disk is viewed 



at intermediate inclination (compare with, e.g.. Fig. 3 in iBeckwith fc Sargentlll993l . which 
illustrates a Keplerian molecular disk model viewed at an inclination of ~45°). The channel 
map closest to the systemic LSR velocity (Kys ~ 2.9 km s~^) indicates a disk position angle 
(PA) on the sky of ~75°, measured east of north. This PA is similar to that inferred for 
the dust disk from the 230 GHz continuum emission detected and resolved by the SMA (see 
§ 13. 3p . We further refine these estimates for disk inclination, PA, and Vgys in § 4.1. 



3.2. CO line profiles 

The SMA ^^C0(2-l) and ^^C0(2-l) emission line profiles, as obtained by spatially 
integrating the line intensity maps within a 6" radius circular aperture, are highly symmetric 
(Figures [2] and |3]). This confirms that t he asymmetric appear ance of the CO line profiles 



obtained with the IRAM 30-m telescope (JKastner et al.ll2008bl ) was a consequence of slight 



mispointing (the 30-m data were obtained at very low source elevation, ~15-20°, due to the 
southerly decination of V4046 Sgr) combined with the spatial resolution of the source by 
the 30-m beam, as opposed to significant density or temperature inhomogeneities intrinsic 
to the molecular disk. 



3.3. 1.3 mm dust continuum emission 

The continuum emission from the dust is partially resolved (Figure H]). Using MIRIAD, 
we fit an elliptical gaussian to the continuum visibilities to estimate the size and shape of 
the emission. The best fit yields major and minor axes of 1.10" and 0.93", respectively, 
with a position angle of 73 ± 5° east of north. This position angle is consistent with that 
inferred from the CO maps (§4.1). There is negligible contribution from the (< 0.1") seeing 
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at 1.3 mm to this Gaussian fit. At a distance of ~73 pc, tlie disk resolved by tlie continuum 
wouid liave a cfiaracteristic radiai extent of about 36 AU (corresponding to tlie FWHM of 
the surface brightness distribution). Both the flux (360 ± 36 mJy) and spatial extent of the 
disk at 1.3 mm are consistent wi th ~40 K blac kbody dust grains located at a distance of 
~40 AU from the central binary ( Kastner et al.ll2008m ). 



4. Modeling 

4.1. CO line maps 

To determine more precise values for the disk parameters, we employed standard mod- 
eling techniques preyiously used to interpret interferometri c maps of molecular line emission 



Qi et al.ll2004l : iRaman et al.ll2006l : iHughes et al.ll2008l). Specifically, w e used the Monte 



(e.g.. 

Carlo radiative transfer code RATRAN (JHogerheijde fc van der TakI 120001 ) to solve for the 
CO level populations at each position within the disk and thereby to generate sky-projected 
emission-line images. These images were then passed to the MIRIAD task uvmodel to sample 
the model images at the same spatial frequencies as our SMA data. A minimum y^ is ob - 
tained by directly comparing the model and observed visibilities (JGuilloteau fc Dutreylll998l ). 
The model geometry is that of a Keplerian disk in hydrostatic equilibrium with truncated 
power laws in density and temperature described (respectively) by n{r) 



nim\ 



lOQAU' 



P and 



Tir) 



Ti, 



100 1 iQOAu ) '^' 'where n(r) refers to the number (mid plane) density of H2 and we 
assume C0:H2 = 10~^. 

There are ten parameters that describe this Keplerian disk model (see Table 2). In 
practice, however, our modeling can constrain only five of these: the central mass, the 
systemic velocity (Vlsr), and the CO disk inchnation (i), PA, and outer truncation radius 
{Rout)- The modeling is relatively insensitive to the other parameters (Tioo, 1, ""-loo; P, and 
inner t runcation radius Rin) because the ^^C0(2-l) emission imaged by the SMA is optically 
thick (JKastner et al.ll2008bl ) and originates from the cooler, outer regions of the disk (the 
^•^00(2-1) emission-line data are too noisy to provide useful constraints on the disk density 
structure parameters nioo and p). Furthermore, there is degeneracy between central mass 
and disk inclin ation; for the forrn e r, we initially adopted the total stellar mass of I.8M0 
determined by IStempels fc Gahml (120041 ) , then varied this parameter so as to explore the 
potential range of i in more detail (§4.2). The values of parameters Tioo, q, "^iioo, P, and 



^The value of ^80 AU quoted in iKastner et al.l ( 2008bl ) is actually a diameter. The semi-major axis is 
-40 AU. 
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Rin (Table 2) were fixed acco rding to results ob t ained from models of the molecular disks 
orbiting other T Tauri stars (JDutrey et al.lll994j . Il998l ). The remaining model parameters 
(Vlsr: P-A., i, Rout) were then varied in turn so as to ascertain best-fit values. The resulting 
set of best-fit model parameters and their formal 1-a uncertainties are listed in Table [2J These 
uncertainties do not account for systematic errors (see below), as they were determined from 
the change in reduced x^ cis the model parameters were varied; a representative contour 
plot of reduced x^ (illustrating its dependence on P.A. and i) is presented in Figure [51 The 
best-fit Keplerian disk model is presented in the form of synthetic velocity-resolved images 
in Fig [6] and as a position-velocity (P-V) diagram alongside the observed P-V diagram in 
Fig El 

The value of Rout = 370 ± 30 AU obtained from the model fitting confirms our ear- 
lier "eyeball" estimate. The systemic velocity we find from fitting the SMA data, 2.92 ± 
0.0 1 km s~^, is consistent with that determined from single-dish, mm-wave CO spectroscopy 
by iKastner et al.l (J2008bl ). The inferred disk inclination of 33.4 ± 0.3° is very similar to the 
value previously assumed in determining the masses of the components of the central binary 
system (see below). The position angle obtained from the best-fit model (76.4 ±0.5°) is con- 
sistent with that derived from the elliptical gaussian fit to the continuum visibility (73 ±5°). 
The truncated power law model generated by these best-fit parameters reproduces well the 
overall emission morphlogies apparent in the velocity channel maps. There is some structure 
in the fit residuals (Fig[6|), however. In particular, the data show an excess of emission 
from the southern portion of the disk relative to the model. The integrated line flux of the 
residuals is ~14% of the ^^CO and model fluxes (see Fig. [2|). Given these systematics, and 
those due to the simplifying assumptions inherent in the molecular disk model, the formal 
errors listed in Table [2l likely somewhat underestimate the uncertainties in the best-fit model 
parameters. 



4.2. Central binary system mass vs. disk inclination 



Estimates of the total mass of the V4046 Sgr close binary, based on radial velocity 
curves obtained from optical spectroscop y of the double-lin ed system, range from 1.55 M© 
(component masses 0.8 6 MfD and 0.69 Mp,; Quast et al.ll2000l) to 1.78 M© (component masses 



0.91 Mq and 0.87 M©; IStempels fc Gahmll2004r). B oth e stimates are based o n an assumed 



binary system inclination of 35° (JQuast et al.l 120001 ). The lStempels fc Gahml (120041 ) compo- 
nent masses appear to be more accurate, given the nearly equal radial velocity amplitudes 
of the two components and the improved agree ment with pre-main sequence evolutionary 
tracks (see disussion in lStempels &: Gahml 120041 ). 



-7- 



To ascertain the constraints placed by the SMA CO maps on both the mass of the 
central binary and the inclination of the disk, we fixed the remaining disk model parameters 
to their best-fit (Table [2]) values, and allowed only central mass and inclination to vary. The 
results for reduced x^, displayed in Fig. El illustrate the inherent degeneracy between mass 
and inclination; lower central masses and higher inclinations yield fits that are comparable 
in quality to those for higher central masses and lower inclinations. Nonetheless, this fitting 
exercise demonstrates that, for a central mass in the range 1.55 to 1.8 M© allowed by the 
optical spectroscopy, the disk inclination is constrained to values between 33° and 37°. Given 
the similarity of these values of disk inclination to that assumed for the central binary in 
determining component masses from stellar radial velocities, it appears that the molecular 
disk is nearly coplanar with the central binary system. Although smaller central binary 
masses (and hence larger disk inclinations) are not precluded by our disk model fitting — 
indeed, we find that x^ reaches a minimum near i ~ 40°, for a central mass of ~1.4 Mq 
— values of central mass < 1.5 M(7^ for this s ystem lead to discord with pre-main sequence 
evolutionary tracks (IStempels fc Gahmll2004l ). 



4.3. Disk gas and dust masses 



To estimate the disk molecular mass from the integrated SMA ^^CO and ^^CO line 
intensities (30.2 Jy km s~^ f or ^^CO and 8 . 7 Jy km s~^ for ^^CO) we adopt the methods 
and assumptions described in iKastner et al.l ( l2008al Jbl) and references therein. In particular, 
we assume optically thin ^^CO emission and C0:H2 and ^^C:^^C ratios of 10~^ and 89, 
respectively; the latter values are highly uncertain. We assume the mean gas temperature is 
similar to that of the dust (37 K; see below). The ratio of the ^^CO and ^'^CO line intensities 
indicates a ^^CO optical depth of ~26. We use this result to infer a dis k molecular mass 
of ~110 Earth masses from the ^^CO line intensity (see Equation 4 in IZuckerman et al. 



20081 ). This gas mass estimate is consistent wit h that obtained for V4046 Sgr on the basis 
of single-dish (IRAM 30-m) data ( Kastner et al.ll2008trl ). We note that the assumption that 
the ^^CO emission is optically thin is supported by the nondetection of C^^0(2-l) emission 
by the SMA, although the large ^^CO optical depth places a lower limit on the ^^CO optical 
depth of ~0.4. 

Adopting a dust temperature of 37 K — obtained from a blackbody fit to the far- 
infrared/submillimeter spectral energy distribution of V4046 Sgr ( jKastner et al.l l2008bl ) — 



Due to a numerical transcription error, the gas mass estimate of ~13 Earth masses stated in lKastner et al 



(J2008bf ) is an order of magnitude too small. 



and a dust opacity of 1.15 cin^ g~^ (extrapolated from an opacity of 1.7 cm^ g^^ at 880 /xm 
assuming a dust opacity slope /?=1; see Zuckerman et al. 2008), we estimate a disk dust 
mass of ~40 Earth masses. This is a factor of ~2 larger than the estimate obtained by 



Kastner et al.l ( l2008bl ): the discrepancy is likely due to the adopted power-law dependence 
of opacity on wavelength. 



5. Summary and Conclusions 

Our ^^C0(2-l) and ^^C0(2-l) emission emission maps of V4046 Sgr reveal a rotating 
disk around this binary system extending out to a radius of ~370 AU. A truncated power- 
law fit to the ^^C0(2-l) visibilities indicate the disk is viewed at an inclination of 33° if the 
total central binary mass is 1.8 M0 (as determined by lStempels fc Gahml 120041 ). However, 
we cannot rule out t he possibility that the central binary mass is as small as 1.55 M© (the 
value determined by lQuast et al.ll2000l ) on the basis of the SMA CO data alone, since disk 
inclinations as large as ~40° (corresponding to a central mass of ~1.4 Mq) are not precluded 
by the modeling. 

While the position angle of the central binary is unknown, the close correspondence 
between the range of disk inclination determined from the modeling (§ 4.2) and the inclina- 
tion of the binary star orbit adopted in previous binary system radial velocity studies (35°; 
Ouast et al.ll2000l ) leads us to conclude that the binary pair and disk are nearly coplanar. The 
disk and orbital planes of circumbinary disk systems are expected to be very similar, because 
the sy stem either formed with such an alignment or subsequently evolved into an aligned 
state ( jMonin et al.ll2006l ). It is therefore significant that the disk and binary inclinations 
indeed appear to be aligned to within a few degrees, in this ~12 Myr-old system. 

We resolved the 1.3 mm continuum emission from V4046 Sgr with the SMA, and 
find that the flux is consistent with 37 K dust located ~40 AU from the central stars 



(JKastner et al.ll2008bl ). We estimate a mass of ~40 Earth masses of dust in the disk, about 
a factor of 3 smaller than the gas mass inferred from the CO line intensities (~110 Earth 
masses, assuming a C0:H2 ratio of 10""^). Hence, either the residual gas mass in the disk 
orbiting this rather "old" binary classical T Tauri system is much smaller than that expected 
for a gas-to-dust ratio of 100, or the disk is severely depleted in gas-phase CO. Indeed, for 
the temperature profile adopted here, T < 20 K for radii ~200 AU, suggesting that much of 
the CO present in the disk may be in the form of icy grain mantles. 



A 



though it appears 



years (JHaisch et al 



most protop l aneta ry disks have dissipated after several million 
2OOII : lUzpenet allboooh . V4046 Sgr, at age -12 Myr, still retains a 
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substantial amount of mass in its circumbinary disk. Indeed, V4046 Sgr is one of only four 
known pre-main sequence stars within 100 pc that possess molecular disks. While V4046 Sgr 
is a close binary system, the other three (TW Hya, HD 141569, and 49 Ceti) do not appear 
to harbor close companions. Previous submillimeter studies have revealed gas-rich disks 



exten ding out to about 200 AU in both TW Hya an d 49 Ceti flQi et al.ll2004J : iHughes et al. 



20081 1 and ~250 AU in HD 141569 flDent et al.ll2005f ). Our SMA imaging demonstrates that 
the V4046 Sgr disk (radius 370 AU) is even larger and, when compared to the disk orbiting 
TW Hya (a single star of slightly later spectral type), also retains a larger gas mass (if we 
apply the same standard assumptions of CO abundance relative to H2 of 10"'^ and ^^C/^^C 
of 89). The twin stars in V4046 Sgr evidently are still actively accreting from this extensive, 
gaseous disk and, likewise, any young Jovian planets within this system may be undergoing 
their final stages of gas accretion. 

Giant planet formation in disks can occur in one of two ways: c ore accretion, with 
timescales of a few xlO^ years (jinaba et al.l l2003l: iLissauer et al.l 120091 ) or disk instability, 
with timescales as short as a few xlO^ years (JBosd I2OOOI ). One expects that long-lived 
disks such as those around TW Hya and V4046 Sgr (ages ~10 Myr) will provide more 
opportunities to form giant planets than shorter-lived disks, regardless of how they form, so 
long as sufficient reservoirs of gas remain present. Meanwhile, studies have shown that disks 
around close (<100 AU) binary system s tend to dissipate faster than those ar ound single 
stars or more widely separated systems ( iBouwman et al.ll2006l : ICieza et al.ll2009[ ). Typically, 
therefore, if giant planets are to form around close binary systems, they must must form more 
rapidly than would be necessary for giant planets forming around single stars. Indeed, there 
is evidence that planets in binary systems with separations <100 AU are biased towards 
higher mas ses, possibly as a natural consequence of rapid planet formation via gravitational 
instability (JDuchenell2010l ). 



V4046 Sgr therefore presents a particularly interesting and important case. Given its 
tight (~9 Rq) separation, the outer portions of its circumbinary disk, where giant planets are 
expected to form, will have been left relatively undisturbed throughout the system's lifetime. 
Only within ~1 AU are dynamical effect s due to the central binary expected to significantly 
affect the disk fjjensen &: Mathieul 119971 estimate that an inner hole out of radius ~0.2 AU 
is necessary to account for the lack of excess infrared emission from V4046 Sgr shortward 
of 10 nm). While about 20% of extrasolar pl anets have been detected in binary systems 
( jRaghavan et al.ll2006l : lEggenberger et al.l 120071 ) . no unambiguous detection of an extrasolar 
planet has been made for a binary with separation <20 AU (though the eclipse timings of 



HW Virginis and CM Draconis suggest that planets may orbit these two binaries; iDeeg et al. 



20081 : iLee et al.ll2009[ ). Hence, there is ample motivation to continue to scrutinize the V4046 
Sgr disk at the highest possible spatial resolution, to see if planets are already orbiting this 
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binary system. 
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Table 1: Observational Parameters 
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Parameter 



i2CO(2~-l) 



13 



C0(2-l) Ci*0(2-1 



Continuum 
l.Smm 



Rest Frequency (GHz) 

Beam Size (FWHM) 

P.A. 

RMS noise (mJy beam~^) 

Peak Flux Density (mJy beam~^) 

Integrated Line Intensity (Jy km s^^) 

Integrated Continuum Flux (mJy) 



230.538 220.399 219.560 225 

2.24" X 1.48" 2.34" x 1.66" 2.35" x 1.67" 1.67" x 1.01' 

go ^o ^o 2° 



87 


76 


76 


0.8 


2640 


790 


<300 


210.5 


30.2 


8.7 


<1 


- 



360 



Table 2: Model Properties 



Parameter 




Value 


Central mass 


(Mo) 


l.S'* 


Inclination (° 


) 


33.4 ±0.3 


Systemic LSR 


'. Velocity (km s"""^) 


2.92 ±0.01 


P.A. (°) 




76.4 ±0.5 


RiniAU) 




4" 


Rout (AU) 




370 ± 30 


riioo (cm-3) 




4 X 10^ '' 


P 




2" 


Tioo (K) 




28^ 


q 




0.6'' 



"See § HJ 

* Adopted values; these are not fit. 

Quoted uncertainties are la. See text for more details. 
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Fig. 1.— SMA velocity-resolved maps of V4046 Sgr in ^^ CO (2-1) (top) and i3cO(2- 
1) (bottom). The contour levels are 0.09 mJy/beamx [3,6,9,...] for ^^C0(2-l) and 
0.08 mJy/beamx[— 2,2,4,6,...] for ^^C0(2-l) with dotted lines representing the negative 
contours. In each frame, the central velocity (in km s~^ with respect t o the Local Standard 
of Rest (LSR)) is indicated at upper left; the position of the continuum source centroid is 
indicated by a small black triangle; and the synthesized beam size, shape, and position angle 
is indicated by the grey shaded ellipse at lower left. 




2 4 

LSR Velocity (km/s) 



Fig. 2. — Line profile of ^^C0(2-l) obtained from the SMA images of the V4046 Sgr disk. 
Overlaid in bold is the best-fit model as described in § HI The dotted line shows the residual 
when the model is subtracted from the data in the visibilities. 
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Fig. 3. — Line profile of CO (2-1) obtained from the SMA images of the V4046 Sgr disk. 
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Fig. 4. — Continuum emission map (left) and deprojected visibility plot (right). The decline 
in the real component of the visibility as a function of baseline is a characteristic of resolved 
structures. The FWHM corresponds to a radial extent of about 36 AU, consistent with the 
characteristic extent of 40 K dust grains in the system. 
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Fig. 5. — Contour maps of reduced x^ for a systemic velocity of 2.92 km s ^ and outer radius 
of 370 AU. The black ellipses indicate 1-, 3-, and 5-0" contours. 
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Fig. 6. — ^^C0(2-l) best-fit model map (top) and residual map (bottom). In each frame, 
the central velocity (km s~^) is indicated at upper left; the position of the continuum 
source centroid is indicated by a small black triangle; and the synthesized beam size 
(2.3" X 1.6" at 8.5°) is indicated by the grey shaded ellipse at lower left. Contours lev- 
els are 0.09 mJy/beamx[— 3,3,6,9,...] for both panels with dotted contours for the negative 
levels. Figure [2] illustrates how the residual flux varies over velocity. The integrated line flux 
of the residuals is ~14% of the ^^CO and model fluxes. 
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Fig. 7. — Position-velocity diagram for the ^^C0(2-l) emission in V4046 Sgr (left) and the 
best-fit Keplerian model (right). Both diagrams were constructed for a P. A. cut of 76.4° 
through the position of the central binary star. 
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Fig. 8. — Contour plot of reduced x^ for disk inclination vs. central (stellar) mass. The black 
contours indicate 1-, 3-, and 5-cr contours; the jagged appearance is due to the coarseness 
of the grid used. The range of central masses covered in the plot is significantly larger 



than that determined vi a studies of stellar radial velocities (1.55-1.8 M©; iQuast et al.l 12000 



Stempels &: Gahml 120041 ) under the assumption of a binary inclination of 35°. 



